Magnetic field induced charge and spin instabilities in cuprate superconductors 
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A d-wave superconductor, subject to strong phase fluctuations, is known to suffer an antiferro- 
magnetic instability closely related to the chiral symmetry breaking in (2+l)-dimensional quantum 
electrodynamics (QED3). Based on this idea we formulate a "QED3 in a box" theory of local in- 
stabilities of a d-wave superconductor in the vicinity of a single pinned vortex undergoing quantum 
fluctuations. As a generic outcome we find an incommensurate 2D spin density wave forming in the 
neighborhood of a vortex with a concomitant "checkerboard" pattern in the local electronic density 
of states, in agreement with recent neutron scattering and tunneling spectroscopy measurements. 



Among the open questions in modern condensed mat- 
ter physics, few have inspired more theoretical effort than 
the emergence of a superconducting state from the doped 
antiferromagnetic (AF) insulator |Q. Recently, using an 
"inverted" approach to the problem [|[ ||, it has been 
shown that AF order arises naturally when the super- 
conducting order in a d-wave superconductor (dSC) is 
destroyed by vortex-antivortex fluctuations (J, H . As we 
shall discuss, the implications of these theories transcend 
the possibility of providing a route to understanding the 
destruction of superconductivity in strongly-underdoped 
cuprates; indeed, they also apply to the problem of local 
field-induced vortices within the superconducting state. 

Recent neutron scattering || ^ and scanning tunnel- 
ing spectroscopy (STS) Q experiments have revealed 
the presence of local AF and charge order in the vicin- 
ity of field- induced vortices. Existing theoretical treat- 
ments §, 0, H|, m of vortex-induced AF ordering 
rely on the proximity of the system to a quantum critical 
point. Within such treatments, it is the suppression of 
the SC order parameter near the vortex cores that leads 
to the nucleation of islands of AF order. Here we present 
an alternative scenario in which the AF order is brought 
about by local quantum fluctuations of a vortex around 
its equilibrium position. In the present theory there is 
no competition between the SC amplitude and AF order: 
the latter arises purely from the presence of vortex fluc- 
tuations and is a genuine low-energy phenomenon taking 
place on lengthscales much longer than the core size. 

It is a well-known fact that the low superfluid den- 
sity in cuprates makes the SC order vulnerable to phase 
fluctuations jl4|, [D| [T(|. This observation has inspired 
theories in which the pseudogap state is modeled as a 
phase-disordered d-wave superconductor |], 17|, [L8|, 
such that the demise of superconductivity is brought 
about by the unbinding and proliferation of the topo- 
logical defects - vortices - in the phase of the SC order 
parameter. It has been pointed out || that fluctuat- 
ing vortices produce a non-trivial Berry-phase interac- 
tion between the quasiparticles of the underlying dSC. 
This interaction is described in terms of a massless non- 
compact gauge field a M , minimally coupled to the Dirac 
fermions representing the low-energy quasiparticle exci- 
tations of the system. Within the theory of Ref. || which 



maps the problem onto (2+l)-dimcnsional quantum elec- 
trodynamics (QED3), it is this interaction that destroys 
the Fermi liquid nature of quasiparticles in the pseudo- 
gap state and ultimately drives the AF instability Q [|. 
Remarkably, both the 'algebraic' Fermi liquid describing 
the symmetric pseudogap phase and the antiferromagnet 
emerge from the same QED 3 theory §J. 

Here, we use the philosophy and formalism developed 
in Refs. || ||, |B| to model quasiparticle excitations in the 
superconducting state in the spatial region close to a sin- 
gle field-induced vortex undergoing fluctuations around 
its equilibrium position. We call this model "QED3 in 
a box" . We note that there exists direct experimental 
evidence that individual vortices indeed undergo signif- 
icant quantum fluctuations p9| . We find that, under 
generic conditions, interactions generated by such fluctu- 
ating vortex lead to local instability of the superconduct- 
ing state which takes form of a 2D incommensurate AF 
spin density wave (SDW) with a wave vector tied to the 
positions of the nodes in the underlying d-wave gap. 

In order to motivate our model for a single vortex we 
first review the treatment of the AF instability in QED3 
and reformulate it in a way that will be more suitable 
for our present purposes. We start from the Euclidean 
QED 3 action S = J cPxCd with 
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describing the low-energy fermionic excitations of a d- 
wave SC coupled to fluctuating vortices represented by 
the gauge field a M ||. Here, ^i(x) is a four component 
Dirac spinor representing the fermionic excitations asso- 
ciated with a pair of antipodal nodes, x = (r, r) denotes 
the space-time coordinate, and 7 M (/1 = 0,1,2) are the 
gamma matrices satisfying {7^,7^} = 25^. The num- 
ber N of fcrmion species is equal to 2 for single-layer 
cuprates; N — 4,6,... for bilaycr, trilayer and multi- 
layer materials. The Lagrangian Cb encodes the dy- 
namics of the gauge field a M and is given by £s[a] = 
^v{q)a^{q)a v {-q) with 



II^(g) = I m a + — -\q\ ) ( 5^ v - 



(2) 



2 



The gauge field mass m a vanishes when vortices are un- 
bound (i.e., in the pseudogap regime or, in the present 
situation, near a single fluctuating vortex) and is finite 
in the superconducting state where vortices appear only 
in tightly bound loops or pairs. 

In the standard treatment Q H the AF order oc- 
curs via the phenomenon of chiral symmetry break- 
ing @ H || in the QED 3 Lagrangian (@). The in- 
stability is signaled by the spontaneous generation of 
fermion mass, ron, which is interpreted in our context 
as the onset of SDW gap [[|, |) for the original Bogoli- 
ubov quasiparticles. The most general, nonperturbative 
treatment of mass- generation in QED 3 obtains as a 
solution of a self-consistent Dyson-Schwinger equation. 
Here we shall follow a slightly simpler route which leads 
to the same result and has the advantage of being more 
easily generalizable to the present problem. In Eq. (|l|) 
we integrate out the gauge field to obtain the following 
fcrmionic effective action: 



(3) 



where Ju{x) = ^f(x)j^(x) is the fermion 3-current and 
D^vix) is the Fourier transform of the gauge boson prop- 
agator D^ v (q) — Henceforth we shall focus 
on a single pair of nodes and thus drop the nodal in- 
dex I. The integrand of the interaction term may be 
written as D^ix — y)Tr[^(y)^(x) r y IJ . i ^(x)^(y)j l/ ] where 
the trace is taken over the spinor indices. This form 
suggests a Hartree-Fock (HF) approach in which we 
decouple the 4- fermion interaction to obtain D^ v (x — 
y)TT[^{y^x)^G Q {x,y) lv ]wit\iG Q {x,y) = (¥(a;)¥(y)) 
and the average is taken with respect to the HF effective 
action to be specified shortly. To make the structure 
of the interaction term more transparent we utilize the 
relative and center of mass coordinates r = x — y and 
R — (x + y)/2 to write it as 



tive action and self-consistency condition: 



S eS = / d A x^{x)~i^id^{x) 



d 3 x / d 3 yJ fli (x)D fil/ (x - y)J v {y), 



d 3 R / d 3 rTv[^(R + )^(R^G (R,r)^] D^(r 



(4) 



where R± = R ± £ . In the uniform system the Green's 
function is independent of R, Go(R,r) = Ga(r). Fur- 
thermore, both Go(r) and D^ v (r) are strongly peaked 
at r — > 0. The dominant contribution to the interaction 
therefore comes from this region and we may write (^) 
as M 



d 3 R^(R)^{R) / d 3 nf ,G {r)-f u D^(r). (5) 



We have dropped the trace since the interaction is pro- 
portional to the unit matrix in the spinor space. 

Inspection of Eq. (||) suggests the following HF effec- 
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d 3 x^(x)(-ff,id^ - im D )^{x), (6a) 



(6b) 



Tr / d'r^GoirhvD^ir). 



The last integral is easily evaluated by going to 
momentum space and Eq. (^) becomes mn = 
(8m£>/iV7r 2 ) ln(A/m£>), where A is the high-momentum 
cutoff. This yields a nontrivial solution 



m D 



= Ae-^ 8 , 



(7) 



in agreement with the classic result of Pisarski . More 
sophisticated treatments pi], based on the Schwinger- 
Dyson equation give a finite critical value of N c above 
which no mass is generated; however, for our purposes 
the level of approximation embodied by Eq. (^) will be 
sufficient. 

We have thus seen that, in a uniform system, fluctu- 
ating vortices lead to the formation of SDW order. The 
challenge we now face is twofold: (i) we must adapt the 
above treatment to the case of a single fluctuating vor- 
tex, and (ii) since we seek to study the commensuration 
effects present in real materials, we must formulate the 
corresponding theory on the lattice. To address (i) let us 
denote by £ v the characteristic length scale over which 
the vortex fluctuates around its classical equilibrium po- 
sition. Within this length scale, the Berry-phase interac- 
tion between quasiparticles (and hence tendency towards 
SDW ordering) will be strong. We model this by taking 
in this region the gauge field to be massless. On the other 
hand at distances well beyond £ v , quasiparticles feel no 
interaction and we model this by gauge field having a 
large mass m a . In particular we take, 



m Q (R) = A ( l ^ 



(8) 



where Ao is an energy scale which we take to be the max- 
imum superconducting gap, |R| is the distance from the 
vortex equilibrium position and n is a positive exponent. 
(We use n = 2 but our numerical calculations below are 
largely insensitive to the exact value of n.) 

To address (ii), (i.e., to put the theory on the lattice) 
we recall that the effective action ([!]) and its HF version 
Eq. (|6a|) descend from a model of a lattice dSC linearized 
near the nodes of the gap. We therefore consider the cor- 
responding lattice Hamiltonian enriched by the "mass" 
term present in Eq. (|6a|), to represent the HF decoupled 
Berry phase interaction. Thus we have, 
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Here $J CT = {c\ a ,c is ), c\ a represents the electron creation 
operator at lattice site i, spin index a = —a, and 



-tij + S i:j (m la - n) 
At- 



Ay \ 

Uj - 5ij(m is - fx) J ' 
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with Uj the tight binding hopping amplitude, Ay the 
SC gap, /i chemical potential, and the local spin 
magnetization representing the mass gap raj) in Eq. (6a). 

We diagonalize Hhf by the generalized Bogoliubov 
transformation c ia = J2 n l u ncr^ihno- + crw* 5 .(r i )7^ (? ], 



where Xno-( r ») = [u n <r(n), av n <r(ri)] satisfy 



(10) 



In terms of the Xna, the self-consistency condition ( pb| ) 
can be written as 



ia = ^af(e n a)V i (rj)ul ia .('Ri+rj)u na (Ri-rj), 
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Vi(r) 



\J dre-^Tr[ ltl D^{T,v) lv ]. (11) 



For m a ^ the above integral cannot be evaluated in 
closed form. However, we find that it can be accurately 
approximated by a simple interpolation formula 



Vi(r) ~ V 



ci 



r(rro a (Rj) + ci)(re 

r. 
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(12) 



where r = |r|, ci = 2/tt and Vo = 16/-/V7T 2 for the case 
of an isotropic Dirac cone (t = A). In the physical case 
t > A the constant Vo will be modified somewhat and 
in what follows we treat it as an adjustable parameter 
of the model measuring the strength of the interaction. 
It is interesting to note that, as seen from Eq. (|l2|), in 
(2+l)D a gauge field mass does not lead to exponentially 
decaying interactions on long length scales. 

To capture the effect of vortex fluctuations on the lo- 
cal superconducting order we solve the eigenproblem fll(i| ) 
numerically on a lattice oi MxM sites and iterate to self- 
consistency using Eq. (|il"|). For simplicity we consider 
only nearest-neighbor hopping amplitudes = t and a 
uniform d-wave gap Ay = ±Ao, with + and — signs 
referring to vertical and horizontal bonds respectively. 
We emphasize that the vortex, fluctuating around the 
equilibrium position at the center of our lattice, enters 
through the position dependent gauge-field mass m Q (R) 
given by Eq. (p|), which in turn enters the potential Vj (r) 
given by Eq. (fL3) . In the spirit of our working philosophy 
that the SDW order arises from the vortex fluctuations 
(and therefore from the gauge field), we neglect at this 
stage any effects of the superflow around the vortex or 
suppression of Ay in the core. Such effects are well un- 
derstood and will be included in a future publication. 
We also neglect the effects of changes in the fermionic 
spectrum due to the onset of SDW on the interaction 
mediated by Berry gauge field Eq. (|l2"l). 

The diagonalizations are performed using standard 
LAPACK routines, which allow us to handle systems 
up to 40 x 40 sites. Typically, 10-15 iterations are 
needed to ensure self-consistency in m ia . We use both 
periodic and free boundary conditions and find that 
they have negligible effect on the results reported be- 
low. Our typical results are summarized in Figs. [I] and 



a) spin magnetization 



b) spin magnetization FT 
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c) staggered magnetization 
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FIG. 1: Magnetization Mi and staggered magnetization Mf 
induced in the vicinity of the fluctuating vortex. Parameters 
used: t = 1, Ao = 0.7, /i = —1.6 resulting in maximum gap 
of 2.3 and average charge density n = 0.62 electrons per site, 
Vo = 1.0 and £ v — 12; periodic boundary conditions. 



^, showing the spatial distributions of the spin magne- 
tization Mi = Xlo- <J ( c L c io')' staggered spin magnetiza- 
tion Mf = Mi(—l) Xi+v \ local electron charge density 

iii = J2cr( c ia Ci c)> anc ^ ener gy integrated local density of 
states (LDOS) S§l(i) = pi{E)dE where Pi(E) is the 
LDOS at site i, as well as their respective Fourier trans- 
forms (FTs). 

Panel (a) in Fig. |l| illustrates the "2D" incommensu- 
rate SDW pattern emerging in the vicinity of a fluctuat- 
ing vortex with an 8 x 8 unit cell containing islands of AF 
order separated by anti-phase domain walls [apparent in 
panel (c)] . The FT displayed in panel (b) reveals that this 
pattern can be thought of as a superposition of four ID 
SDWs with wave vectors Qsdw = tt(1±<5sdw> 1±(5sdw)> 
Ssdw = \- The size of <5sdw is doping dependent: 
it shrinks with increasing fi and vanishes at half filling 
(/i = 1), giving rise to perfectly commensurate AF SDW. 
We also find that for fi < —1.7 the SDW becomes very 
weak for reasonable values of coupling Vo: overdoped 
samples are less susceptible to AF instability. 

According to general symmetry arguments, a spatial 
modulation in the spin density generates a modulation 
in the charge density, Srii oc Mf. For our 2D SDW 
pattern this implies that the corresponding CDW will 
have a unit cell with half the area, rotated by 45° rela- 
tive to the unit cell of Mf. Indeed, panels (a) and (b) 
of Fig. U confirm this general expectation, showing a 
"checkerboard" CDW at principal wavevectors Qcdw = 
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a) charge density 



b) charge density FT 





c) integrated LDOS 
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d) integrated LDOS FT 
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FIG. 2: Local charge density n; and integrated LDOS So.'lCO 
for the same parameters as Fig. [I]. 



7r(±5cDW,0),7r(0, ±<5 C dw) with <5 C dw = \- A similar 
checkerboard pattern arises in the integrated LDOS and 
is displayed in panels (c) and (d). 

Our findings of a checkerboard pattern in LDOS are 
consistent with the recent STS experiments performed 
on Bi 2 Sr 2 CaCu20 8+ a (BSCCO) crystals §. Our pre- 
diction is that the period of the pattern should increase 
with underdoping and that the effect should vanish in 
the overdoped samples. Also, if the observed LDOS pat- 
tern is associated with electron density modulation in 
a single CuO layer, we predict that the corresponding 
neutron scattering peaks should be found at wavectors 



Qsdw = tt(1 ±| ( 1± t). We note that neutron experi- 
ments H 0| on La 184 Sr .i6CuO4 (LSCO) show peaks at 
different fc-space positions, 7r(l ± \, 1) and 7r(l, 1 ± j). 
Although the findings of STS and neutron experiments 
are generally cited as being mutually consistent our anal- 
ysis above indicates that this is not necessarily so: for a 
genuine 2D SDW illustrated in Fig. [j] determination of 
the corresponding CDW must consider the interference 
terms which cause the apparent 45° rotation of the latter. 
In the absence of neutron measurements on BSCCO we 
see two possible resolutions of this difficulty. First, it may 
be that the CDW pattern is truly 2D and neutron scat- 
tering on BSCCO would find a pattern illustrated in Fig. 
|l|(b). Second, it could be that STS sees an incoherent 
superposition of two orthogonal ID CDWs originating 
in two CuO layers comprising the BSCCO bilayer. This 
would explain the x — y anisotropy reported in Ref. § 
and the neutron pattern would be consistent with that 
observed in LSCO. Within our simple model such ID so- 
lutions have slightly higher energy than the 2D solutions 
reported above but it is possible that in a more complete 
model (using e.g. a more realistic band structure) the 
situation will be reversed. 

To summarize, we have presented a "QED3 in a box" 
theory for field induced spin and charge instabilities in 
cuprates, driven by local fluctuations of a pinned vor- 
tex. Without any need to fine-tune parameters we find 
LDOS patterns in detailed agreement with the tunneling 
data on BSCCO § and we relate them in a plausible 
way to existing neutron experiments |], 0. Other mod- 
els [|[ 10, [ll], [l2|, [l3| rely on the local suppression of 
the superconducting order parameter and would there- 
fore predict similar instabilities in the vicinity of impuri- 
ties, grain boundaries, and sample edges where thus far 
no such effects have been observed. 
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